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Abstract. Selected aspects of the description of neutron-induced fission in ^*'^Pu in the 
framework of the nuclear energy density functional theory at finite temperature are presented. 
In particular, we discuss aspects pertaining to the choice of thermodynamic state variables, the 
evolution of fission barriers as function of the incident neutron energy, and the temperatures of 
the fission fragments. 



1. Introduction 

Nuclear fission remains one of tlie most complex physics phenomena in nature. Powerful 
phenomenological models, finely tuned to a wealth of precise experimental data, have given us a 
good qualitative, and at times quantitative, understanding of the phenomenon. However, these 
models of fission lack a firm connection with the theory of the nuclear force. This jeopardizes 
their ability to supply the high-precision predictions currently needed to understand the later 
stages of the formation of elements in supernova, or to simulate new prototypes of nuclear 
reactors. 

Some of the basic concepts of a microscopic theory of fission were already laid out in 
the eighties by the theory group at CEA Bruyeres-le-Chatel in France [H [2]. Starting 
from an effective interaction between nucleons that embeds (in a phenomenological way) in- 
medium many-body correlations, the nuclear self-consistent mean- field theory provides a generic 
framework to construct ever more accurate representations of the nuclear wavefunction [3j. To 
describe fission, a number of collective variables must be introduced. The potential energy 
surface in the collective space thus defined serves as a basis to perform a time evolution of the 
nuclear wave packet in the framework of the time-dependent generator coordinate method. As 
evidenced from the few applications reported, the method is very promising |4j. 

However, at the time of its inception, the necessary computing power was lacking so that 
this approach could not yet compete with more empirical models. Today, the fast development 
of leadership class computers and interdisciplinary collaborations such as the NUCLEI project 
open new perspective In addition, recent progress in deriving quality energy functionals for 
fission applications [6] and of understanding the quantum mechanics of scission suggests the 
microscopic approach to fission has enough potential to become predictive. 

In this work, we discuss some aspects of the microscopic theory of fission in the poster-child 
example of the neutron-induced fission of ^^^Pu at low energies. After a brief reminder of the 
theoretical framework, we show the evolution of the fission barriers as function of the nuclear 
temperature, and discuss some of the practical aspects related to such simulations. 



2. Theoretical framework 

Our formalism is based on the nuclear energy density functional theory (DFT) with local energy 
functionals of the Skyrme type [8] for the particle-hole channel, and contact density-dependent 
interactions with mixed volume and surface character for the particle-particle channel [9j. The 
nucleus is described at the Hartree-Fock-Bogoliubov approximation: its ground state is a quasi- 
particle vacuum that depends on a number of collective variables q = {gi}i=i,...,Ar, such as, e.g. 
the expectation value of multipole moments or angular momentum. The set of the A'^ collective 
variables defines a point q in the collective space. For each point in the collective space, the 
requirement that the energy be minimal with respect to variations of the generalized density 
leads to the HFB equations; solving the latter defines the one-body density matrix and the 
pairing tensor, and thereby all observables of interest. 

The finite-temperature extension of the HFB formalism (FT-HFB) was introduced in nuclear 
physics more than 30 years ago |10j . We only recall that in the finite-temperature HFB theory, 
the ground-state is a statistical superposition of pure quantum states, characterized by a density 
operator D. Adopting the HFB approximation for the density operator leads to the finite- 
temperature HFB equations: they are formally equivalent to their form at zero-temperature, 
only the expression of the density matrix and pairing tensor is different. Most importantly, the 
Wick theorem still applies, meaning that all observables are computed from the trace (in the 
given representation) of the corresponding operator and the one-body density matrix. 

In addition to the density operator, macroscopic thermodynamical concepts are also used to 
describe the system at finite temperature. In particular, several thermodynamical potentials are 
available to describe the nucleus, each coming with its set of state variables. The two most useful 
are (i) the internal energy Eiy^S^X), where V is the (constant) nuclear volume, S the entropy 
and X any additional extensive state variable needed to characterize the nucleus (expectation 
value of multipole moment, for example) and (ii) the Helmholtz free energy F(y, T, X) with 
F = E — TS and T the temperature. The FT-HFB equations are obtained by minimization 
of the grand potential at temperature T, and naturally yield the free energy in the isothermal 
representation (constant T). Passage to an isentropic description (constant S) is sometimes 
needed, as discussed below. 

3. Induced fission and finite-temperature formalism 

In the description of induced fission in the finite-temperature DFT framework, the first step is 
to compute the potential energy surface of the nucleus as a function of the chosen collective 
coordinates and the nuclear temperature. In the following, we work in a four-dimensional 
collective space characterized by the expectation values of the axial, Q20-, and triaxial, Q221 
quadrupole moment, mass octupole moment Q30 and hexadecapole moment Q40. All calculations 
have been performed with the SkM* parametrization of the Skyrme functional pTj. In 
the particle-particle channel, the proton and neutron pairing strengths were fitted to the 
experimental 3-point formula for the odd-even mass differences in ^^Op^. a cutoff of E^nt = 60 
MeV limits the number of quasi-particles taken into account in the calculation of the densities. 

Calculations were performed with the DFT solvers HFBTHO [l2] and HFODD [13j. In both 
codes the solutions to the HFB equations are expanded in the one-center harmonic oscillator 
(HO) basis. In HFBTHO axial and time-reversal symmetry are assumed so that solutions are 
labeled by the projection of the angular momentum on the z-axis. By contrast, HFODD is fully 
symmetry-unrestricted, and can, in particular, describe triaxial and parity-breaking shapes. The 
two programs have been benchmarked against one another and agree within a few eV for an 
axial configuration |T2]. In practice, the HFB wavefunctions were expanded onto A'states = 1100 
HO basis states with contributions from up to A^max = 30 oscillator shells. 



3.1. Evolution of fission harriers 

In figure [U we show the evolution of the free energy F = E — TS as function of the expectation 
value Q2Q of the axial quadrupole moment. All three other degrees of freedom are locally 
minimized: in practice, the first barrier has a non-zero value of Q22 while the octupole moment 
is non-zero past the second barrier. Note the gradual disappearance of the barriers as the 
temperatures increases. 




Figure 1. Evolution of fission barriers in ''Pu as a function of temperature along the most 
probable fission path. Plain lines curves with a symbol show the free energy at constant 
temperature (isothermal process), the dashed- line curves next to them show the corresponding 
internal energy at constant entropy (isentropic process). 



In applications of induced fission, the primary motivation has to do with the properties of the 
fission fragments rather than those of the fissioning compound nucleus. Of particular importance 
are the charge and mass distributions of the fragments, their total kinetic energy (TKE), and 
their excitation energy. The latter is a major input in the reaction codes that model the neutron 
and gamma spectrum generated during fission. Estimates of TKE and excitation energy require, 
in principle, to work in the isentropic representation. The Maxwell relations of thermodynamics 
state that, for any (extensive) collective coordinate qi, 



dF 
dqi 



dE 

dqi 



Therefore, variations of internal energy with respect to deformation at constant entropy are 
equal to the variations of free energy at constant temperature }14j. 

Prom the set of curves {F{q,Tk)}k obtained directly from the FT-HFB calculations, one 
can, by numerical interpolation, reconstruct the E{q., S) for any value S of the entropy. The 
procedure is based on a standard spline interpolation. Finite-temperature HFB calculations 



produce at each point q in the cohective space the quantities E(T)\q and S{T)\q, where from 
we easily obtain E{S)\q for any desired value of S. When normalized to the same deformation 
point, the curves E{q)\s and F{q)\T should be strictly equivalent within the interpolation errors. 
This is verified in figure [U where all dashed lines are obtained for the values of entropy at the 
top of the first barrier. 

3.2. Temperature of the fragments 

Until now, macroscopic-microscopic approaches assumed either the same temperature in the 
two fragments |15] . or a ratio of temperatures proportional to the ratio of the fragment masses 
T1/T2 ~ A1/A2 |16] . The validity of the Maxwell thermodynamical relations in a practical 
case suggests a method to predict microscopically the nuclear temperature of each individual 
fragment after scission. 
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Figure 2. Temperature along the most probable fission path of "Pu as extracted from the 
entropy as function of the internal energy. 



As seen in the previous section, the isothermal and isentropic representations are equivalent. 
The advantage of the latter is that it involves only extensive variables, and is, therefore, more 
amenable to the description of the splitting of the compound nucleus into two subsystems. The 
scission of the compound nucleus, characterized with an internal energy E and an entropy S, 
yields two fragments, with respective internal energy and entropy {Ei,Si) and (£'2,5'2). Since 
the entropy is an extensive variable, we have S = Si + S2 (here, we have to assume implicitly 
that the interaction energy between the two fragments is 'small enough'). The next step is to 
take advantage of the thermodynamical relation 



For each fragment n = 1,2, we can perform a separate FT-HFB calculation constrained on the 
deformations q^""^ of the fragment at scission. This gives us the functions 5^") Taking 
the derivative with respect to E^"^^ gives the temperature T^"-* of the fragment n. While the 
value of both the internal energy and entropy in this separate FT-HFB calculation should be 
the same as the ones extracted from the compound nucleus (again: if the interaction energy can 
be neglected), the temperature may be different, as it is related to the derivative dS/dE. 

The method just outlined will be applied to estimating the temperatures of the fragments in 
the induced fission of ■^^'^Pu in a forthcoming publication. As a preliminary step, we estimate in 
figure 2 its numerical accuracy by comparing the exact value of the temperature, as set in the 
FT-HFB calculation, with the numerical estimates obtained from the S{E) function. Overall, 
the accuracy is of the order of 50 keV at worse, owing to interpolation errors. 

4. Conclusions 

Recent advances in the nuclear energy density functional theory, combined with the constant 
increase of computing power, have enabled significant progress toward the development of a 
microscopic theory of nuclear fission. In this short article, we have presented early results 
related to the evolution of fission barriers with nuclear temperature and the equivalence between 
isentropic and isothermal representations. 
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